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改进 PSO 算法 在 多 无 人 机 协同 任务 分 配 中 的 应 用 


葛 ” 硕 ， 喜 小 平 
(中 国 矿 业 大 学 ,江苏 徐州 221116) 


摘 要 : 针对 多 无 人 机 协同 任务 分 配 越 来 越 复杂 的 问题 ， 采 用 一 种 改进 的 阶层 分 级 粒子 群 优 化 算法 (HGIWPSO) 获 得 最 
优 分 配方 案 。 首先， 根据 粒子 适应 度 值 将 种 群 动 态 划 分 为 不 同 三 个 阶层 , 依据 不 同 阶层 粒子 特性 选择 合适 的 学 习 模 型 ， 


并 引入 独立 权重 思想 调节 惯性 权重 大 小 ， 平 衡 算法 全 局 与 局 部 搜索 能 力 ， 提 高 算法 性 能 ; 然后 ， 建 立 协同 多 任务 分 配 
问题 模型 ， 采 用 多 余 负 载 竞 拍 方案 减少 非法 劣 解 ， 通 ww 解决 实际 
分 配 问题 。 实 验 结果 表明 ， 该 算法 能 够 有 效 解 决 复杂 约束 条 件 下 多 无 人 机 协同 任务 分 配 问题 ， 得 到 最 优 分 配 序列 ， 具 


有 一 定 的 理论 以 及 实际 意义 。 
关键 词 : 无 人 机 ; 阶层 分 级 ; 任务 分 配 ; 粒子 群 算法 ; 实数 编码 
中 图 分 类 号 : TP301.6 doi' 10.3969/j.issn.1001-3695.2018.05.0303 


Application of Improved PSO algorithm in multi UAV cooperative task allocation 


Jiang Shuo, Yuan Xiaoping 
(China University of Mining & Technology, Xuzhou Jiangsu 221116, China) 


Abstract: In order to solve the increasingly complex problem of multi-UAV collaborative task allocation, used an Improved 
hierarchical particle Swarm optimization algorithm (HGIWPSO) to obtain the optimal allocation scheme. Firstly, according to 
particle fitness, the population dynamics were divided into three different classes, and selected the appropriate learning model 
according to the characteristics of different classes of particles, and introduced the independent weight idea to adjust the 
weight of the inertia weight that improved the performance of the algorithm. Then, established the model of collaborative 
multi task assignment problem, and adopted the redundant load auction Scheme to reduce illegal solutions. Established the 
mapping relationship between the particle and the actual allocation Scheme through real number coding to solve the actual 
allocation problem. Finally, the experimental results show that the algorithm can effectively solve the problem of multi-UAV 
coordinated task assignment under complex constraints and obtain the optimal distribution sequence. It has certain theoretical 
and practical significance. 
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乘客 , 由 此 将 问题 简化 来 解决 UAV 航路 规划 任务 问题 ,文献 [2] 


中 将 UAV 任务 分 配 问题 分 解 转换 为 经 


的 MTSP 和 NFO 问题 
多 UAV 协同 决策 分 配 问 题 规模 大 小 和 作战 过 程 中 将 不 同 求解。 文献 3，4] 对 多 个 UAV 建立 约束 条 件 以 及 目标 函数 来 建 


任务 进行 耦合 求解 等 因素 决定 了 其 复杂 性 。 目 前 国内 外 对 UAV ” 立 MILP 模型 求解 ,文献 [5] 分 析 各 种 类 型 任务 之 间 的 约束 关系 ， 


任务 分 配 问题 的 研究 主要 集中 数学 建 模 和 分 配 模型 的 求解 优化 。 将 UVA 的 总 航程 和 所 需要 的 时 间 进 行 优化 , 建立 CMTAP 模型 
这 两 个 方面 。 现 阶段 建 模 方面 的 主要 成 果 包 括 了 车 辆 路 径 问 题 。 求解 。 
模型 (vehicle routing problem, VRP)、 多 旅行 商 问 题 模型 (multiple 对 模型 求解 算法 主要 有 禁忌 搜索 、 遗传 算 法 千 吕 下 开 到 引用 源 ,- 鲁 如 | 


traveling salesman problem，MTSP)、 混 合 整数 线性 规划 模型 ”到 用 汤 。 (genetic algorithm，GA)、 蚁 群 算法 千 ! 不 并 到 引用 浙 。(ant 
(mix-integer linear programming，MILP)、 网 络 流 优化 模型 colony optimization, ACO) 和 粒子 群 优化 算法 等 丰 找到 引用 源 , 错 洪 ! 未 找到 相 


(network flow optimization, NFO) 和 多 UAV 


(multiple processors resource allocation, CMTAP) 等 。 文献 [1] 
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主要 研究 方向 为 电能 质量 检测 与 控制 . 物 联网 技术 . 


办 同 任务 分 配 模 型 用 浙 。 (particle swarm optimization, PSO) 等 。 文 献 [12] 设 置 禁忌 


将 ”来 表征 已 有 局 部 最 优 位 置 ， 在 求解 过 程 中 不 断 更 新 与 测试 局 部 
UAV 看 做 VRP 模型 中 的 车 辆 , 将 侦察 目标 看 做 VRP 模型 中 的 。 最 优 值 ， 最 终 求 得 最 优 解 。 文 献 [13] 设 计 遗 传 算法 对 多 UAV 


表 
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同 分 配 问 题 进行 求解 ， 并 与 分 支 界 定 法 和 随机 搜索 方法 进行 比 
较 表 明 其 具有 较 好 的 时 间 效 率 。 文 献 [14] 通 过 二 进 制 矩阵 编码 
实现 UAV 的 离散 化 编码 ， 然 后 以 交叉 变异 策略 加 强 粒 子 多 样 
性 , 将 PSO 算法 有 效 的 运用 于 求解 多 UAV 侦察 任务 分 配 问 题 。 
文献 [15] 提 出 一 种 多 UAV 任务 分 配 与 航 迹 规划 整体 控制 架构 ， 
采用 与 协同 系统 相 匹配 的 粒子 改进 PSO 算法 解决 任务 分 配 问 


六 


沼 
> 


上 述 求解 方法 都 能 在 一 定 程度 上 解决 UAV 协同 任务 分 配 
问题 ,但 是 都 存在 求解 算法 复杂 或 者 求解 效率 较 低 等 问题 ,PSO 
算法 参数 较 少 、 原 理 简单 、 搜 索 速 度 较 快 ,适合 用 于 求解 UAV 
协同 任务 分 配 问 题 。 本 文 提 出 的 阶层 分 级 策略 ， 针 对 不 同 阶层 
粒子 的 特性 给 出 相应 更 新 方式 ， 增 强 种 群 多 样 性 ， 提 高 算法 的 
收敛 精度 以 及 搜索 效率 ， 并 对 粒子 采用 实数 编码 方式 简化 算法 
求解 UAV 协同 任务 分 配 问题 ， 取 得 较 好 的 效果 。 


1 ”模型 建立 


1.1 模型 描述 

本 文 涉 及 的 主要 环境 为 模拟 二 维 战场 环境 ， 假 设 战 场 中 的 
目标 任务 集合 为 了 = 他,7…Ty} ，N, 为 需要 执行 的 目标 任务 总 
数 。 对 于 vTieT ,我 们 通过 一 个 六 元 组 来 标志 任务 目标 的 属性 ， 
<Tid7,Pos7 ,Stater ,Value’ ,D’,Thr” > 为 相应 的 六 元 组 。 昌 中 每 个 任 
务 的 标记 都 是 唯一 确定 的 ， 用 Tid 标志 。 在 二 维 战场 空间 中 用 
Pos" 来 标志 其 坐标 位 置 ， srarer 代表 当前 任务 目标 的 状态 ， 
Value” 代表 UAV 攻击 指定 的 任务 目标 价值 量 , 也 就 是 目标 任务 
的 价值 ，D" 标志 的 是 敌 方 目标 的 威胁 范围 区 域 ，7Thr" 代表 的 是 
UAYV 攻击 敌 方 目标 对 UAV 自身 的 威胁 系数 , 也 就 是 UAV 执行 
该 任务 被 击落 的 概率 ， 

本 文中 UAV 的 集合 为 VY={V,Vy…Ww,} ，N 代表 的 是 UAV 
的 数目 。UAV 属性 用 个 七 元 组 表示 ， 
< Uid", Pos" ,State” ,Value" ,At,Ts,L> 为 相应 的 七 元 组 。 其 所 代表 的 
内 容 和 目标 任务 类 似 ， 分 别 为 执行 目标 的 UAV 编号 、UAYV 在 
二 维 战场 上 的 坐标 位 置 、 当 前 的 状态 、UAV 本 身 拥有 的 价值 、 
执行 目标 UAV 的 攻击 能 力 、 任 务 集合 以 及 最 大 的 任务 执行 载 
和 荷 。4rs[00 ,也 即 表 示 UAV 攻击 任务 目标 的 时 候 能 够 击毁 任务 
目标 的 概率 ， 工 表示 的 是 UAV 能 够 执行 的 最 大 任务 数量 ， 其 
中 1<|75|<L。 
1.2 代价 函数 

UAYV 对 敌 方 目标 进行 打击 的 过 程 中 , 已 方 被 敌对 目标 击毁 
的 可 能 尽 可 能 小 。 设 第 i 架 UAV 执行 完 任务 j 后 的 存活 概率 为 
Pp; ， 则 p=1-Thn 。 因 此 ， 一 架 UAV 完成 N 个 任务 所 需 付出 的 
代价 值 如 式 (1) 所 示 。 


Thr” el[0,1) 。 


Ni 
Ci=Value’ *(1-[ {0-7hn)) (1) 


/j=l 


其 中 : Value? 为 UAV i 的 价值 ，Thr 为 UAV 被 敌对 目标 /7 击 
的 概率 。 
在 UAV 与 敌对 攻 生 


落 


二 < 


LI 


目标 都 处 于 较为 理想 的 状态 下 时 ,己方 
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执行 攻击 任务 的 UAV 距离 目标 任务 点 的 距离 大 小 ,决定 了 UAV 
航程 代价 的 大 小 ，UAV 相对 越 接近 目标 其 航程 代价 将 会 越 小 ， 
因此 目标 被 距离 较 近 的 UAV 执行 的 几率 相对 更 大 。 航 程 代价 
的 表达 式 如 (2) 所 示 。 


C2=dy/ dax (2) 
其 中 4; 为 UAV i 与 目标 j 的 欧式 距离 ，4d wm 表示 所 有 UAV 
相距 目标 j 的 最 远 距离 ， 即 d wx = maxi(di) 。 
攻击 收益 指 的 是 UAV 执行 完成 攻击 任务 , 击毁 任务 目标 获 
得 目标 估 值 ， 该 项 评价 标志 在 对 目标 进行 优化 决策 中 引导 最 终 
结果 趋向 于 作战 能 效 最 大 化 ， 如 式 (3) 所 示 。 
C3=Value, * 4 (3) 
基于 上 文中 的 三 项 评价 标准 并 非 统 一 的 度量 单位 ， 具 有 不 
可 公 度 性 而 无 法 进行 简单 的 累加 综合 计算 ， 需 要 通过 一 些 标 准 
化 变换 对 相应 的 评价 标准 进行 处 理 。 本 文选 取 线 性 尺度 变换 法 
对 UAV 价值 和 目标 估 值 对 比 效益 指标 实现 规范 化 操作 ， 如 式 
(4)(5) 所 示 。 


_ Value’ 


max Value” 
ieN, 


*a-] [0-7)) 由 


Value' 


3 二 后 
max Value, 
jeN, 


A (5) 


多 UAV 协同 分 配 实质 上 是 一 个 优化 多 目标 的 问题 , 本 文通 
过 引用 和 借鉴 线性 加 权 和 法 将 复杂 的 优化 多 目标 的 问题 转化 为 
相对 简单 的 优化 单 目标 的 问题 。 通 过 设置 相应 的 权重 向 量 来 平 
衡 各 个 因素 对 分 配 结果 的 影响 程度 。 其 中 权重 向 量 表示 为 


SS 


且 满足 ws(0D ，> ,wmw=1。 多 UAV 任务 分 


四 =(@Ob 2 03) ， 


配 问题 转换 为 求解 式 (6) 所 示 最 优 问 题 ) 9， 
minf -Seertoc-ocyn) (6) 


其 中 : 决策 变量 轧 为 0-1 变量 ， 思 =1 表示 UAV i 执行 任务 j 。 


约束 条 件 如 式 (7)~ 式 (9) 所 示 。 
Dw <l, VvieJ (7) 
MN, 
Xi <L viel (8) 
AL 
YY < No Amin(N,, SL) (9) 


在 上 面 所 述 的 相关 约束 条 件 中 ， 式 (07) 表示 每 一 项 指定 的 攻 
击 任务 只 能 由 一 架 UAV 去 完成 ; 式 (8) 指 的 是 每 架 UAV 所 能 执 
行 的 任务 载荷 限制 的 大 小 ， 比 如 一 个 UAV 的 任务 载荷 为 3， 那 
么 该 架 UAV 最 多 允许 执行 3 个 任务 ; 式 (9) 表 示 的 是 整个 机 群 
能 够 执行 任务 的 总 量 大 小 ， 也 就 是 总 共 能 够 执行 的 任务 数量 的 
大 小 ， 通 常 这 个 约束 条 件 都 能 够 满足 ， 可 以 保证 每 一 个 任务 都 
能 够 被 执行 ;UAV 与 目标 的 下 标 集 合 定义 分 别 为 7 人 由 2…X,) 和 
JA{l,2,..., N,} 。 
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2 ”改进 的 PSO 算法 


2.1 基本 PSO 算法 
基本 PSO 算法 的 数学 描述 如 下 和 针 ! 末 找到 引用 浙 。 错误 ! 未 找到 引用 源 。， 
粒子 群 的 初始 种 群 规模 定 为 M ， 粒 子 维 数 为 上 维 ， 在 :时 
刻 粒 子 i 在 解 空间 的 坐标 定义 为 ,Xi =Gi,xo… 巩 ) ， i= 上 2,.…,M ， 
其 速度 大 小 决定 了 粒子 每 次 在 解 空间 飞行 距离 大 小 ， 用 
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粒子 在 搜索 过 程 中 对 于 自身 经 验 和 群体 经 验 都 需要 借鉴 ,所 以 ， 
对 于 中 层 阶 级 粒子 的 速度 更 新 方式 保持 标准 PSO 算法 的 速 
更 新 公式 不 变 如 式 (14) 所 示 。 

W =v +on(py -XW +on(g; —X) (14) 
那么 表明 粒子 为 下 层 阶 级 粒子 。 
对 于 下 层 阶级 粒子 ， 其 搜索 位 置 已 经 远 远 偏离 了 最 优 位 置 ， 自 
身 搜索 经 验 参考 价值 不 大 ， 所 以 对 于 下 层 阶级 粒子 速度 更 新 方 


赣 


QC) 如 果 3M /4<NX, < NX ， 


VW =04, 咏 ,… 苏 ) 表示 ， 则 粒子 1 在 时 刻 1 的 第 jw =12.… 四 维 空 


间 中 的 速度 以 及 位 置 更 新 公式 如 式 (10)~(12) 所 示 。 

Vi = Ov + cn(py 一 2 ) 十 cj7(8 1 一 2 ) (10) 
Vy Vi > 
让 -| Tk (GD 
辣 
t _ ji-l ft 

XH 二 XN 十 六 (12) 
其 中 : @ 为 惯性 权 值 (inertia weight)，c 和 < 为 加 速 因子 


(acceleration constants )， 志 和 是 在 [0,1] 内 的 两 个 随机 数 ， 粒 
子 在 解 空间 进行 搜索 的 过 程 中 ， 粒 子 的 速度 会 被 限定 在 
大 的 搜索 范围 之 内 ， 速 度 最 大 值 用 vi 表示， 限定 粒子 在 搜索 
过 程 中 的 范围 ， 改 善 算法 性 能 。8; 表示 人 进 代 过 程 中 群体 的 历史 
最 优 位 置 ， 即 全 局 极 值 gr ， 而 P; 表示 的 则 是 粒子 自身 搜索 过 
程 中 的 最 优 位 置 ， 即 个 体 极 值 Pw 。 

2.2 阶层 分 级 的 PSO 算法 

PSO 优化 过 程 中 ， 当 代 全 体 粒 子 的 适应 度 值 大 小 用 集合 
了 ={ 有 记 …f} 表示 , 其 中 严 代 表 粒 子 的 总 数 , 代表 的 是 第 i 个 
粒子 的 适应 度 值 大 小 。 根 据 计算 所 得 适应 度 值 对 全 体 粒 子 进行 
升序 排列 (此 处 适应 度 值 越 小 代表 粒子 位 置 越 优 ), 形成 新 的 粒 
子 序列 集合 =f, 为,…,%} ,通过 判断 粒子 距离 第 一 个 粒子 距离 
的 远近 将 种 群 划分 为 不 同 的 三 个 阶层 。 设 定 全 局 最 优 粒子 重新 
排序 后 的 序列 号 为 NX。， 具 体 的 阶层 分 级 策略 如 下 详 述 : 

a) 如 果 和 X。 SNX; <M 14， 那么 表明 粒子 为 上 级 阶层 粒子 。 
其 中 ，ANX; 表示 重新 排序 后 第 i 个 粒子 的 序列 号 ，M 代表 总 的 
粒子 数量 ， 通 过 对 粒子 重新 排序 ， 粒 子 的 序列 号 标志 了 其 距离 
当前 最 优 粒 子 的 远近 。 对 上 级 阶层 粒子 ， 其 自身 搜索 已 经 处 于 
最 优 值 附 近 ， 随 着 算法 不 断 迭 代 ， 粒 子 自身 越 来 越 靠 近 最 优 位 
EE ， 这 个 时 候 整 个 群体 的 经 验 对 于 自身 来 说 已 经 不 具备 太 大 的 
学 习作 用 ， 所 以 粒子 在 迭代 更 新 过 程 中 忽略 群体 学 习 经 验 ， 主 
要 考虑 自身 学 习 经 验 模 型 。 对 此 ， 我 们 对 于 上 级 阶层 粒子 在 原 
有 速度 更 新 公式 上 提出 新 的 速度 更 新 公式 ， 适 应 粒子 自身 的 特 
点 ， 达 到 更 好 的 搜索 寻 优 的 目的 ， 上 级 阶层 的 粒子 速度 更 新 公 
式 如 式 (13) 所 示 。 


个 最 


多 =OV tan(py -区 ) (13) 

其 中 ，1t 代 表 粒 子 的 当前 迭代 次 数 。 
b) 如 果 必 14<NX; <3M /4， 那 么 表明 粒子 为 中 层 阶 级 粒子 。 
对 于 中 层 阶 级 的 粒子 ， 其 自身 的 搜索 位 置 处 于 优 务 位 置 之 间 ， 


式 ， 采 用 群体 学 习 模型 ， 其 速度 更 新 如 式 (15) 所 示 。 
= 十 CR(8) 一 地 ) (15) 
同时 ， 因 为 下 层 阶 级 粒子 已 经 远 远 偏离 了 最 优 解 的 位 置 ， 
为 了 能 够 让 该 层级 的 粒子 能 够 和 优秀 粒子 进行 经 验 信息 交流 ， 
让 下 层 阶级 的 粒子 和 全 局 最 优 粒 子 以 一 定 的 概率 ? 进行 交叉 操 
作 ， 促 进 粒 子 的 信息 交流 。 随 机 选择 两 个 点 42 ， 满 足 
a,be[l,Dl,a<b， 其 中 DD 表示 粒子 维 数 , 交叉 二 者 之 间 的 维 数 形 
成 新 的 粒子 。 

PSO 算法 在 迭代 后 期 容易 失去 多 样 性 ,算法 陷入 局 部 最 优 ， 
为 解决 该 问题 ， 算 法 中 引入 随机 扰动 速度 跳出 局 部 陷阱 。 定 义 
第 i 个 粒子 到 当前 最 佳 粒子 的 维度 距离 为 


(16) 


此 


式 (17) 表 示 维 度 距离 的 标准 


w=- (17) 
其 中 : 4 表示 维度 距离 的 期 望 值 。 
当 粒 子 的 状态 满足 VV< 只 时 (其 中 办 为 定义 的 标 称 值 )， 判 
定 算法 在 迭代 过 程 中 己 经 陷入 局 部 最 优 陷阱 。 此 时 给 与 一 个 随 
机 速度 扰动 ， 帮 助 粒子 摆脱 局 部 最 优 ， 继 续 进行 迭代 搜索 ， 随 
机 扰动 速度 更 新 公式 如 式 〈18) 所 示 。 


Ov + on(py -x )+V *rand(), NX, < NX < 

1 | 让 M 3M 

六 二 Ov ton(py -Wh ton(g) -+ rand0, F< NX < (18) 
Ov +cn(8; 一 大 ) + trandO, <NX, < NXy 


结合 本 文 改进 策略 的 特点 ， 选 取 独 立 惯性 权重 调节 粒子 全 
局 与 局 部 寻 优 能 力 ， 如 式 (19) 所 示 乞 下 并 到 引用 浙 


和 
@ = on -our = Oo ) EP (19) 


4 为 指数 变化 


Ba 


其 中 @w 为 初始 权重 值 ， 
子 ， 计 算 方式 如 下 所 示 : 


%w 为 结束 权重 值 ; 


万 Go) 
fit,, (20) 


arctan(log 


KE =$+n* 
A 


其 中 pn 为 固定 常量 ， 应 C6) fi 分别 代表 粒子 i 的 适应 度 值 和 
平均 适应 度 值 。 通 过 对 数据 进行 数学 变换 取 值 是 为 了 数据 尽量 
满足 假设 ， 并 且 让 波动 相对 稳定 。 
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站 六 到 引用 浙 ,。 实验 通过 统计 适应 度 误 差 值 的 均值 与 标准 差 以 及 而 
究 适 应 度 值 随 迭 代 次 数 的 变化 曲线 来 评价 算法 改进 后 的 优化 性 
3.1 算法 性 能 能 ， 取 500 次 的 平均 值 ， 仿 真实 验 结果 如 下 所 示 。 
使 用 测试 函数 来 测试 算法 性 能 , 标准 测试 函数 如 表 1 所 示 乱 
表 1 实验 采用 的 测试 函数 及 其 相关 属性 


3 HGIWPSO 算法 求解 任务 分 配 问 题 


十 


名 称 表达 式 初始 化 范围 
Sphere f1(x)= > 2 [-100,50 
il 
Rosenbrock f2(x) = (000 一 + 1) [-2.048,2.048] 
i=] 
1 电 1 
Ackley f3(x) =e+20—20exp( > exp(5 2 cos(277)) [-32.768,16] 
Griewank f 4(%) >” x Leox Fn #1 [-600,200] 
Rastrigin f5(x) = > (x -10cos(2zx)+10) [-5.12,2] 


i=l 


表 2 算法 对 比 数据 (10 维 ) 


测试 函数 PSO SAPSO CLSPSO HGIWPSO 
而 mean 0.0098 4.3384e-14 7.6352e-59 6.1328e-99 
std 0.0278 3.0633e-13 5.3983e-58 4.0220e-98 
区 mean 7.1355 0.6437 0.7973 0.5267 
std 2.0210 1.4863 1.6105 1.3520 
二 mean 2.3077 0.7969 1.5510 0.6543 
std 0.9419 0.8525 0.8652 0.8252 
而 mean 0.1649 0.1404 0.1603 0.1376 
std 0.1090 0.0709 0.0876 0.0833 
mean 12.4199 8.8153 16.3571 6.7856 
std 4.8183 4.5882 7.0616 2.8898 
表 3 算法 对 比 数据 (30 维 ) 
测试 函数 PSO SAPSO CLSPSO HGIWPSO 
- mean 4.7752e+02 2.4555 0.0338 0.0118 
std 1.8189e+02 1.2729 0.4943 0.1705 
区 mean 1.5298e+02 29.8197 12.3050 3.1052 
std 59.5632 9.4281 16.8979 4.2535 
芝 mean 8.8051 3.6526 6.5942 2.8480 
std 1.1737 0.8119 1.3359 0.6245 
有 mean 1.3858 0.0926 0.0389 0.0144 
std 0.1493 0.0434 0.0285 0.0141 
mean 89.2660 71.4109 96.5783 54.3095 
std 16.3387 20.2922 25.0982 15.3889 
| | 10 9 
10 
雪 划 
En 二 二 加 0 
一 人 一 HGPSO 10 
-po0| 一 HGIWPSO 
10 
0 1000 2000 3000 4000 5000 0 1000 2000 3000 4000 5000 
迭代 次 数 和 迭代 次 数 


对 


图 1 Sphere 函数 收敛 对 比 曲 线 (10 维 ) 


2 Sphere 函数 收敛 对 比 曲 线 (30 维 ) 
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0 1000 2000 3000 4000 S5000 0 1000 2000 3000 4000 5000 
迁 代 次 数 迁 代 次 数 
图 3 Rosenbrock 函数 收敛 对 比 曲线 (10 维 ) 图 4 Rosenbrock 函数 收敛 对 比 曲线 (30 维 ) 


0 1000 2000 3000 4000 5000 0 1000 2000 3000 4000 5000 
迁 代 次 数 迁 代 次 数 
图 5 Ackley 函数 收敛 对 比 曲线 (10 维 ) 图 6 Ackley 函数 收敛 对 比 曲线 (30 维 ) 


PSO 
IWPSO 
HGPSO 
HGIWPSO 


0 1000 2000 3000 4000 5000 0 1000 2000 3000 4000 S5000 
迭代 次 数 迭代 次 数 
图 7 Griewank 函数 收敛 对 比 曲 线 (10 维 ) 图 8 Griewank 函数 收敛 对 比 曲线 (30 维 ) 
10 3 
六 
E10 
咖 
10" 
0 1000 2000 3000 4000 5000 0 1000 2000 3000 4000 S5000 
迭 代 次 数 迭代 次 数 
图 9 Rastrigin 函数 收敛 对 比 曲线 (10 维 ) 图 10 Rastrigin 函数 收敛 对 比 曲线 (30 维 ) 
图 表 数 据 可 以 很 明显 地 看 出 ， 不 管 是 对 于 高 维 情况 还 是 数量 等 于 粒子 维 数 , 粒子 的 整数 部 分 用 来 标志 执行 任务 的 UAV 
低 维 的 情况 ，HGIWPSO 算法 相 比 于 其 他 对 比 改 进 算法 能 够 跳 编号 ， 相 同 的 整数 部 分 则 代表 着 一 架 UAV 执行 任务 的 编码 序 
出 局 部 最 优 位 置 ， 取 得 更 好 的 优化 精度 ， 较 好 的 收敛 速度 。 误 列 ， 小 数 部 分 以 升序 排列 决定 UAV 执行 攻击 任务 的 先后 次 序 。 
差 和 均 方 差 均 优 于 其 它 的 改进 算法 ， 相 比 其 他 算法 具有 更 好 的 同时 引入 拍卖 机 制 优化 减少 非法 务 解 ， 保 证 负载 约束 条 件 能 


性 能 ， 证 明 改 进 的 可 行 性 和 有 效 性 。 
3.2 HGIWPSO 算法 实现 任务 分 配 
本 文采 用 实数 编码 方式 进行 编码 。 


满足 。 
加 入 拍卖 机 制 的 任 
人 初始 化 粒子 群 ; 


务 分 配 算法 执行 流程 如 下 : 


设 定 被 攻击 任务 目标 的 
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b) 对 初始 化 粒子 进行 解码 ， 如 果 满 足 约束 条 件 则 转 到 下 一 ” 配 中 ， 更 具有 优势 ， 能 够 跳出 局 部 位 置 ， 最 终 寻 找到 最 优 位 置 

步 ， 如 果 不 满足 则 按照 拍卖 算法 的 任务 协调 分 配方 法 进行 调整 (2.7473, 1.8132, 1.4449, 3.6773, 2.7441, 1.7111, 3.8870) 

修正 ， 重 新 编码 粒子 ， 更 新 粒子 位 置 ; 让 执行 任务 的 总 代价 最 小 为 0.8433， 得 出 最 优 分 配方 案 如 表 6 
c) 计 算 各 个 粒子 的 适应 度 值 ; 所 示 。 并 且 在 170 代 左 右 就 成 功 收敛 到 最 优 位 置 ， 收 敛 速 度 较 
d) 更 新 粒子 速度 和 位 置 ; 快 。 证 明 本 文 提出 的 改进 算法 应 用 于 多 UAV 协同 任务 分 配 上 
©) 判 断 收敛 准则 是 否 成 立 ， 如 果 成 立 ， 给 出 分 配方 案 ,， 算  ” 面 效果 较 好 ， 改 进 策略 可 行 ， 具有 一 定理 论 和 实际 意义 。 


法 结束 ， 如 果 不 成 立 转 入 b) 继 续 执 行 。 
4 ”结束 语 
给 出 三 架 执 行 任务 的 无 人 机 以 及 需要 被 执行 攻击 任务 的 目 
标 参 数 如 表 4 和 5 所 示 生 下 到 引用 泊 。 氏 并 1 术 坟 到 引用 浙 。。 本 文 针 对 多 UAV 协同 任务 分 配 复杂 多 耦合 的 特性 ， 给 出 
攻 4 目标 参数 可 行 建 模 方案 ， 简 化 问题 进行 求解 。 采 用 阶层 分 级 策略 并 结合 
目标 编号 位 置 价值 杀伤 概率 独立 权重 的 改进 PSO 算法 很 好 的 解决 了 任务 分 配 问 题 。 改进 算 
Tl (55, 50) 20 0.50 法 采用 阶层 分 级 和 独立 权重 思想 ， 不 同 于 现 有 的 改进 策略 ， 大 
T2 (60，70) 20 0.55 多 针对 整个 群体 ， 而 是 通过 阶层 分 级 充分 体现 不 同 阶层 粒子 的 
T3 (65，30) 40 0.60 寺 性 , 并 由 独立 权重 策略 平衡 个 体 粒子 全 局 以 及 局 部 搜索 能 力 ， 
T4 (75，20) 100 0.85 能 够 很 好 的 增强 粒子 的 多 样 性 ， 让 粒子 能 够 更 加 容易 跳出 局 部 
T5 (80，50) 40 0.70 最 优 ， 加 快 搜索 速率 ， 很 好 的 弥补 了 算法 的 缺陷 ， 并 相对 于 现 
T6 (78, 70) 50 0.75 阶段 相关 改进 算法 效果 更 为 显著 。 将 改进 算法 用 于 求解 多 UAV 
T7 (75, 90) 30 0.80 协同 任务 分 配 问题 ， 采 用 实数 编码 方式 降低 问题 的 求解 难度 ， 
表 5 无 人 机 参数 并 通过 负载 竞拍 方式 解决 实际 负载 约束 问题 , 取得 较 好 的 效果 。 
编号 位 置 价值 任务 载荷 ” 杀伤 概率 
UAV1 (15，45) 50 3 0.60 SN 
UVA2 (20，50) 60 4 0.70 [1] Zorlu O. Routing unmanned aerial vehicles as adapting to capacitated 
UAYV3 (15, 55) 40 2 0.80 vehicle routing problem with genetic algorithms [Cl]// Proc of International 
为 验证 算法 在 任务 分 配 应 用 中 的 有 效 性 ， 进 行 仿真 对 比 实 Conference on Recent Advances in Space Technologies. 2015: 675-679. 
验 。 算 法 进行 仿真 实验 时 参数 如 下 : 种 群 大 小 为 40， 粒 子 维度 [2] Km M H, Baik H, Lee S. Response threshold model based uav search 
为 7 维 ， 最 大 迭代 次 数 300， 战 术 权 重 分 别 取 (0.2.0.4,.0.4)， 国 planning and task allocation [J]. Journal of Intelligent & Robotic Systems， 
定常 量 分 别 取 0.73 和 1.02， 惯 性 权重 最 大 值 和 最 小 值 分 别 为 2014, 75 (3-4): 625-640. 
0.9 和 0.4。 取 50 次 结果 的 平均 值 ， 仿 真实 验 结果 如 下 所 示 。 [3] 叶 媛 媛 ， 闵 春平 ， 朱 华 勇 ， 等 . 基于 整数 规划 的 多 UCAYV 任务 分 配 问题 
表 6 任务 分 配 结果 研究 [加 . 信息 与 控制 ，2005，34 (5): 548-552. (Ye Yuanyuan, Min 
无 人 机 任务 分 配 计划 Chunping, Zhu Huayong, etal. Research on multi UCAYV task allocation 
UAV1 T3_T6_T2 based on integer programming [J]. Information and Control, 2005, 34 (5): 
UAV2 TS 一 Tl 548-552. ) 
UAV3 T4_T7 [4] 张 庆 杰 . 基于 一 致 性 理论 的 多 UAV 分 布 式 协同 控制 与 状态 估计 方法 
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